Radiationless Travelling Waves In Saturable Nonlinear Schrodinger Lattices 
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The longstanding problem of moving discrete solitary waves in nonlinear Schrodinger lattices is 
revisited. The context is photorefractive crystal lattices with saturable nonlinearity whose grand- 
canonical energy barrier vanishes for isolated coupling strength values. Genuinely localised travelling 
waves are computed as a function of the system parameters for the first time. The relevant solutions 
exist only for finite velocities. 
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Recently, the topic of discrete solitons (intrinsic lo- 
calised modes) in photorefractive materials has received 
significant attention, see for a review. This interest 
was initialized by the experimental realization of two- 
dimensional periodic lattices in photorefractive crystals 
in which solitons were observed 0, Q . Further work 
has revealed a wealth of additional coherent structures 
such as dipoles, quadrupoles, soliton trains, vector, neck- 
lace and ring solitons, see e.g. 0,0,0, El- Since photore- 
fractive materials feature the so-called saturable nonlin- 
earity, these results for periodic lattices have spawned a 
parallel interest in genuinely discrete saturable nonlinear 
lattices 0, 0, 0| . One particularly intriguing result of 
theses studies is that the stability properties of the lo- 
calized modes are substantially different from their reg- 
ular discrete nonlinear Schrodinger (DNLS) analogs. In 
DNLS it is well-known jl^J that site-centered lo- 

calised modes are always stable, while intersite-centered 
modes are unstable (and are only stabilized in the con- 
tinuum limit where these two branches degenerate into 
the well-known continuum sech-soliton of the integrable 
cubic NLS). For the photorefractive nonlinearity (i.e., 
the so-called Vinetskii-Kukhtarev model originally pro- 
posed in 15j and revisited in 0, 0, depending 
on the coupling strength, the inter-site centered modes 
may have lower energy than their onsite counterparts. 
Hence, one should expect that the ensuing sign rever- 
sal of the so-called Peierls-Nabarro (FN) energy barrier 
A£' = Eis — Eos (where the subscripts denote inter- 
site and onsite respectively) should cause an exchange 
between the linear stability properties of the two modes. 

A related, even more fundamental, question in non- 
linear lattice models of DNLS type is whether exponen- 
tially localized self-supporting excitations that move with 
a constant wave speed can exist, so called moving discrete 
solitons. For continuum models this question is in some 
sense trivial since the equations posed in a moving frame 
remain of the same fundamental type. Yet for lattices, 
passing to a moving frame leads to so-called advance- 
delay equations that are notoriously hard to analyse. One 
recent (negative) result in this direction for the so- 



called Salerno model, shows that, starting from the inte- 
grable Ablowitz-Ladik equation, mobile discrete solitary 
waves acquire non-vanishing tails as soon as parameters 
deviate from the integrable limit. Hence exponentially 
localised fundamental (single humped) moving discrete 
solitons cannot be constructed. While these results set- 
tle a long-standing controversy, see e.g. 0,0,0, they 
are also somewhat unsatisfactory since they do not give 
conditions under which moving discrete solitons might 
exist for generic, non-integrable lattices. 

In the arguably simpler problem of kinks in Frenkel- 
Kontorova lattices (and, more generally, in discrete 
Klein-Gordon models), genuinely travelling fundamental 
('charge one') topological solitons are known not to oc- 
cur unless there is a "competing" nonlinearity that causes 
vanishings of the FN barrier; see e.g. |2fl l2lj. Drawing 
an analogy from this, for (the more complicated, com- 
plex field) DNLS models with pure cubic nonlinearity, 
the FN barrier never disappears, so one should not ex- 
pect genuinely localised moving solitons (in some sense 
anticipating the above-mentioned negative result). How- 
ever the present saturable model has a nonlinear term in 
its denominator, which, under Taylor expansion, is effec- 
tively cubic for low intensities but exhibits saturation for 
large intensities. This effective "competition" between 
nonlinear behaviors leads to the possibility of the vanish- 
ing of the FN barrier, and hence suggests this model as a 
good candidate for finding moving discrete solitons. Such 
a search is the theme of this Letter. We note that such a 
connection was also partially discussed in [i^, based on 
the properties of the stationary localized modes. 

We start by examining the existence of stationary 
localized modes as the strength of inter-site coupling 
is varied, initializing our computations from the anti- 
continuum limit, where explicit solutions are available. 
A linear stability analysis of the solutions is carried out, 
and we find that the regime of stability-instability al- 
ternation does not coincide with the vanishings of the 
FN barrier but rather with the vanishings of the grand- 
canonical free-energy of the system G = E — AP, where A 
is the frequency (chemical potential) of the solution and 



2 



P is its P norm. After identifying this sequence of trans- 
parent points, we proceed to obtain, for the first time in 
discrete nonhnear Schrodinger type models, genuinely lo- 
calized, single-humped, uniformly traveling solutions in 
their vicinity. This is done b y ex amining travelling so- 
lutions as embedded solitons [23, such that while reso- 
nant with the linear spectrum in the travelling frame, the 
tail amplitude can vanish for isolated, definitely non-zero 
values of the speed, for a given set of system parameters. 
This is reminiscent of recent work in the Klein-Gordon 
context [2^ , where the so-called Stokes constant is shown 
to vanish for isolated values of the speed, indicating the 
presence of travelling localised solutions. 

The model of interest, as discussed above, will be a 
discrete lattice with photorefractive nonlinearity, namely, 

TO„ = -£A2U„ + — — y T^Un. (1) 

Here, e denotes the coupling strength between adjacent 
sites and /3 is the coefficient of the nonlinearity (propor- 
tional to the voltage in photorefractive crystals) and can 
be scaled out for our purposes; A2W„ = m„+i-|-m„_i — 2u„ 
denotes the discrete Laplacian. We will therefore set 
/3 = 1 in what follows. We seek stationary solutions of 
the form u„ — e~*^*w„, starting from the so-called anti- 
continuum (AC) limit e = 0. In that limit the principal 
solutions of interest arc the onsite solution with w„„ = 
±y/ (1/A) — 1, (all other sites bearing vanishing ampli- 
tudes) and the intersite-centered solution with m„q — 
Uno+i = ±-\/(l/A) — 1 and zero amplitudes elsewhere. 
We initialize these exact profiles from the AC limit and 
subsequently use continuation in e to solve the station- 
ary equations for the wave profile via Newton's method 
for finite e. We also examine the linear stability of the 
obtained modes using linearization around a solution pro- 
file v° of the form: u„ = e"*^* [v° + S (a„e^* + 6„e^**)] , 
and resolving the ensuing matrix eigenvalue problem en- 
suing to O(^) for the eigenvalue A and the eigenvector 
pair (a„,6fj (where * stands for complex conjugation). 
In what follows, we will also use the (Hamiltonian) en- 
ergy of solutions, 

E^J2 - ""I' + + ' (2) 

n 

and the Z^-norm, P — J2n l^nP- 

Our stationary calculations as a function of the inter- 
site coupling (rather than of the frequency as in 0) are 
shown in Fig.^ Panels A and B show the key imaginary 
eigenvalues of the problem, while panel C shows the max- 
imal real eigenvalues. A solid line is used in the figures to 
represent the inter-site centered (IS) mode, and a dashed 
line the onsite-centered (OS) mode. Observe that for 
small values of e the OS mode is stable while the inter- 
site mode is unstable. Consider the OS mode. At e « 0.1, 
an eigenvalue emerges from the band edge of the contin- 
uous spectrum at A = zA, and approaches the origin of 
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FIG. 1: (Color Online) Panels A and B show the key imagi- 
nary eigenvalues for the onsite (dashed line) and the intersite 
(solid line) mode as a function of e (A on a linear scale, B 
on an exponential scale). The band edge of the continuous 
spectrum is at A = 0.5. From there bifurcates an eigenmode 
for £ > 0.1 (for the onsite case) which arrives at the origin 
of the spectral plane for e — 0.25 and becomes real. On 
the contrary for the intersite solution a previously real eigen- 
mode exits as imaginary for e > 0.25. These modes alternate 
again for e = 0.445, e = 0.7 etc. The real parts of the corre- 
sponding eigenvalues are shown in the panel C. Panel D shows 
log(|Ai5|) between onsite and intersite modes, and by dashed 
line the quantity log(|AG|), where G = E — \P (see text for 
details) . 



the spectral plane. This feature also occurs in the pure 
cubic DNLS but the eigenvalue only arrives at the origin 
of the spectral plane as e — > oo (the continuum limit). 
Here the eigenvalue arrives at the origin for e ~ 0.25, be- 
comes real for 0.25 < £ < 0.445, then becomes imaginary 
again for 0.445 < e < 0.7, then real again, etc. Further- 
more, this alternation of stability occurs hand-in-hand 
with the alternation of stability of the IS mode, which 
starts out unstable, becomes stable, then unstable again, 
with the transitions occurring at precisely the same s- 
values. Note from the figure that stability alternation 
occurs, not at the points of vanishing of the FN energy 
barrier; but when the grand- canonical energy barrier of 
the model AG = A{E - AP) between the IS and OS 
modes vanishes. This provides an alternative viewpoint 
to the results of , which is also rather natural from the 
point of view of statistical mechanics for a system with a 
conserved "particle number" {P norm); see e.g. |2^. At 
these very points, the relevant eigenvalue (pertaining to 
translation) crosses zero, instantaneously restoring an ef- 
fective translational invariance in the model. Thus these 
are transparent points where there is a one-dimensional 
family of localised solutions, obtained from each other by 
translation, of which the OS and IS modes are just two 
examples. 

Intuitively, one might expect that the transparent 
points represent parameter values at which genuinely 
travelling localised solutions of the photorefractive lat- 
tice might bifurcate with zero wave speed. Hence we 
shall look for solutions to Eq. ^ near such points using 
a travelling wave substitution Un{t) = il)ij>' ~ ct)e'*^'^"~^*-' . 
This leads to the advance-delay equation in the travelling 
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the tail of the sohtary wave. At best this could happen in 
a codimension-one set in parameter space, that is along 
discrete curves lying within the single resonance bands of 
Fig. 121 Solutions of this type are sought using a pseudo- 
spectral method, see 0,1131 ■ To this end, we set 



N 

tp{z) — flj cos(cjjz) 



ibj sm{LL!jz) 



(5) 



FIG. 2: Intersections of the single and double root conditions 
(see text) for A — 0.5 and varying c and e. The shaded area 
shows where there is more than one branch of linear waves. 
Values indicate the number of roots of Q. Subplots show 
left and right hand sides of I^J for {e,c) — (1,3), one root 
(i), (3,3), three roots (ii) and (3.5,0.6), seven roots (iii). (iv) 
displays the overlapping of the bands for small e and c, only 
the first six bands have been shown. 



co-ordinate z — n ~ ct 

-icij'{z) = (2e - A)tp{z) - e{ij{z + l)e''' + 'ipiz - l)e-''') 
1 + \i;{z)\ 

where ' denotes differentiation with respect to z. Equa- 
tion iPJ is rotationally invariant, therefore the transfor- 
mation ip{z) = ■0(z)e~*'^^ can be used to obtain the same 
equation as above for ijj and A = fcc + A. Therefore, 
k = will be used hereafter in © without loss of gener- 
ality. The linear spectrum is obtained by substitution of 



into the linear part of Eq. Q 
cA + A - 1 = 4£sin^(-^) 



(4) 



It is interesting to note that the asymptotic behaviour 
of the two sides of Eq. for small A is such that a 
resonance with the linear spectral bands is unavoidable. 
However, in the hope of obtaining a localized solution 
such resonances need to be minimized and hence we need 
to investigate parameter regions where Eq. has only 
one root, A > 0. This can be efficiently done |2l| by 
considering the conditions for double roots, namely c = 
2esin(A), which allow us to distinguish among regions 
with different root multiplicity. The results are plotted in 
Fig. [3 from which it can be seen that for each e > 0, the 
number of spectral bands increases to infinity as c ^ 0. 
This result serves to counter the naive expectation that 
near the transparency points, travelling solutions of small 
speed may exist. 

The best hope for finding genuinely localised solutions 
is then within parameter regions of Fig.|21in which there 
is a single resonance with linear waves. Such solutions 
would represent embedded soliton structures |2ll . |23| for 
which the radiation mode component exactly vanishes in 



with u!j = to transform |j3J) into a system of alge- 
braic equations in the long finite interval [—L/2, L/2]. 
aj,bj G R are the coefficient of the Fourier series, for 
which the ensuing algebraic equations are solved at the 
collocation points Zm = 2{n+i) ■ "^^^ solution is obtained 
by means of the Powell hybrid method 26] , with an error 
tolerance of 10""'^'^. From the resulting aj,bj's, ip{z) can 
be reconstructed and numerically continued using AUTO 
[irf with an error tolerance of 10~^° to investigate the ef- 
fect of varying parameters on the solution shape and its 
tail amplitude. Generally, the thus obtained solutions 
will be weakly nonlocal ones with non-zero oscillatory 
tails. However, appending an additional tail condition 
for the solution, yielding a signed measure A of the am- 
plitude of the tail, we can obtain truly localised solutions 
as isolated zeros of this function A. Given our choice of 
phase in the substitution ||SJ) that the real component of 
Tp is odd about z — L/2 (and vice versa for the imagi- 
nary component), a good choice for such a tail amplitude 
measure is A = Im(-0(^)). 

Continuation of weakly nonlocal solitary wave solu- 
tions as e varies is shown as a function of 1/e in Fig.OKtop 
left). Note the topological distinction between branches 
the 'S '-shaped branches that contain true zeros of A and 
the 'U' or 'n'-shaped branches that do not [l^. Once 
the isolated points where the tail disappears are found, 
these can be continued in any pair of desired parame- 
ters, e.g. {e, c) by appending the condition A = to the 
numerical problem. Thus, we can construct the paramet- 
ric diagram of the existence of such genuinely travelling 
(with constant prescribed speed), yet truly localized (i.e., 
radiationless due to A = 0) solutions. A regular sequence 
of branches so-obtained are found as e increases, the first 
three of which are shown in the top right panel of Fig.|3 
Notice that each branch terminates (as it should) at the 
edges of the multi-spectral bands. Nevertheless, at the 
low-c end the branches are close to straight lines which 
if continued to zero wave speed would hit the transpar- 
ent points. Thus, while there is no actual bifurcation as 
such of intrinsic localised modes from these transparent 
points, their presence has a strong influence. 

The middle two panels of Fig. 13 show solution profiles 
on the first (lowest e) branch of localised solutions. Note 
how the profile delocalizes when the edge of the multi- 
resonance band is reached. Any attempt to continue zero- 
tail solutions into the band resulted in lack of convergence 
for large enough L and N. The lower two panels show 
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FIG. 3; (Top left) Continuation of weakly nonlocal solitary 
waves for various values of 1/e showing three zeros in A for 
c = 0.7, A = 0.5, L = 60. Zeros of A at e « 0.76, 1.02, 1.36. 
(Top right). Continuation of the 3 zeros of A, varying e and 
c with A = 0.5. Circles on the c = axis indicate the trans- 
parency points. Middle panels show the continuation of the 
(real and imaginary parts of the) solution of branch I for val- 
ues of the speed. Examples of direct integration of the so- 
lution of branch I with e = 0.911396, c = 0.894153 (bottom 
panels) . The space time contour plot of the solution modulus 
(bottom left panel) , and the modulus before (solid) and after 
(dashed) 100 time steps (bottom right panel) are shown. 



the result of direct integration of a representative solution 
on branch I, clearly illustrating their genuinely travelling 
nature, and indicating that they are stable against the 
perturbation introduced by numerical discretization. 

In conclusion, to our knowledge, this is the first time 
genuinely localised travelling waves have been found in a 
discrete lattice of nonlinear Schrodinger type. While it is 
true that our construction has relied on numerical contin- 
uation, we have demonstrated the robustness of the solu- 
tions we have obtained, and have provided a new ratio- 
nale as to why they should exist. This rationale is based 
on the existence of transparent points for stationary lo- 
calised modes when a generalised Peierls-Nabarro barrier 
vanishes. However, we have also shown that this is not 
enough; one also needs the wave speed to be sufficiently 
large to overcome the multi-spectral bands. Within each 
band, exponentially localised solutions occur as embed- 
ded solitons which correspond to regular zeros of an ap- 
propriately defined measure of the tail amplitude. 



We anticipate that the principal features of the model 
problem studied (competing nonlinear terms for the ex- 
istence of transparent points and sufficiently high wave 
speeds to overcome the multi-spectral bands) should be 
generic features a variety of similar Hamiltonian lattices. 
Indeed, we have found similar results for a pure cubic- 
quintic nonlinearity. More realistic models of photore- 
fractive crystals should be expected to behave similarly, 
giving the possibility of experimental observation of the 
waves we have constructed. There, the robust, radiation- 
less nature of the new moving pulses may be particularly 
appealing. Another relevant question is the extension of 
the present ideas to two-dimensional settings jl^ , where 
deciding which directions of wave propagation within the 
lattice can lead to radiationless pulses will also be of in- 
terest. 
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